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Abstract 

Nonlinear plane acoustic waves propagating through a fluid are studied using 
Burgers' equation with finite viscosity. The evolution of a simple N-pulse with 
regular and random initial amplitude and of pulses with monochromatic and noise 
carrier is considered. In the latter case the initial pulses are characterized by two 
length scales. The length scale of the modulation function is much greater than the 
period or the length scale of the carrier. With increasing time the initial pulses are 
deformed and shocks appear. The finite viscosity leads to a finite shock width, which 
does not depend on the fine structure of the initial pulse and is fully determined by 
the shock position in the zero viscosity limit. The other effect of nonzero viscosity 
is the shift of the shock position from the position at zero viscosity. This shift, as 
well as the linear time, at which the nonlinear stage of evolution changes to the 
linear stage, depends on the fine structure of the initial pulse. It is also shown that 
the nonlinearity of the medium leads to generation of a nonzero mean field from an 
initial random field with zero mean value. The relative fluctuation of the field is 
investigated both at the nonlinear and the linear stage. 
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1 Introduction 



The propagation of finite amplitude sound waves is of fundamental interest in nonlinear 
acoustics. In the simplest model of propagation in fluids these waves are described by the 
well-known Burgers' equation (plane waves) |I],3 or modifications of Burgers' equation, 
which are called generalized Burgers' equations (cylindrical and spherical waves) [@-[5|]. 
In studies of nonlinear wave propagation an important problem is to find the waveform 
of the asymptotic wave at long time after the preparation of the initial wave or at long 
distance from the source emitting the wave. In the first case the asymptotic wave is 
called the old-age wave and is an asymptotic solution of an initial value problem of some 
wave equation and in the latter case the asymptotic wave is a solution of a boundary 
wave problem. The asymptotic wave is damped both by absorption and by shock wave 
dissipation and the asymptotic wave is therefore, for plane waves, described by the linear 
diffusion equation. Because the linear diffusion equation describes the attenuation of high 
frequency waves, the asymptotic behaviour is determined by the spectrum of the waves of 
low frequencies. This spectrum is the result of the distortion of the wave at the nonlinear 
stage. 

The decay of a wave to the old-age waveform is very different for periodic signals 
and pulse perturbations. For the periodic signal the old-age waveform is an exponen- 
tially decaying harmonic wave. The pulse perturbation has continuous spectrum and the 
amplitude of the pulse decays according to a power law. 

For sinusoidal and N-wave initial perturbations the old-age solutions have been studied 
in several papers for plane, cylindrical and spherical waves with use of both analytical 
and numerical methods [p|-|r5|]. The dependence of the amplitude constant of the old-age 
waveform on the parameters of the initial perturbations has been found. 

The aim of the present paper is to investigate the asymptotic behaviour of complex 
pulses, such as modulated or random waves. Plane waves are studied, which means that 
we will use the original form of Burgers' equation 0. Burgers' equation has an exact 
solution, which is found by reducing it to the linear diffusion equation by mean of the 
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Hopf-Cole transformation fig, |17[. Because of the existence of an exact solution it is 
relatively easy to find the old-age waveform developing from simple initial perturbations 
like periodic signals and N-waves in the plane wave case. However, for signals with complex 



structure it is far from trivial to find their evolution even for plane waves |18-f20|. For 



a regular signal with fractal structure [|T^,|20[ an unusual sequence of stages of evolution 
may appear: the nonlinear stage may succeed the linear stage. At the nonlinear stage 
this wave may decrease more slowly than a periodic wave or a simple pulse. 

Still more complicated behaviour is shown by a random signal. The solution of Burgers' 
equation with random initial conditions is often called Burgers turbulence. Numerous 



papers are devoted to investigations of this problem |22|-p7|. In the case of vanishing 
viscosity the continuous random initial field is transformed into a sequence of straight lines 
with some slope and with random locations of the shocks separating them. Due to the 
merging of the shocks the internal scale of the turbulence increases and the random wave 
decreases more slowly than the periodic signal. The decay rate depends on the behaviour 
of the initial energy spectrum of low frequency [pTjpol and is also sensitive to the statistics 
of the potential of the initial velocity field [0,|26|,|27|. The asymptotic behaviour of the 
field in the case of finite viscosity also depends strongly on the statistical properties of the 
initial field. In the case where there are large scale components in the initial spectrum the 
nonlinear stage never transforms to the linear regime of evolution |21 |. In the opposite 
case the final evolution depends on the tail of the initial probability distribution |23| . 

In the present paper we consider the evolution of complex pulses which are character- 
ized by two scales: Iq - the inner scale of the carrier, Lq - the scale of the modulation, 
and the condition Lq » Iq. For such signals the generation of a low- frequency com- 
ponent or a non-zero mean field takes place. The case of vanishing viscosity has been 
investigated [p9|, pT|1 . It has been shown that, after multiple merging of the shocks, the 
inner structure disappears and finally a finite pulse ends up with an N-wave. This N-wave 
is fully described by the positions of its zero and shocks, which are determined by the 
potential ipo( x ) of the initial velocity v (x): v' (x) = —ip (x). It has also been shown 
that, for a pulse with random carrier, the fluctuations of the positions of the zero and the 
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shocks are relatively small for Iq « Lq. 

In the present paper we will consider the influence of finite viscosity on the nonlinear 
stage of the evolution of the pulse and investigate its old-age behaviour. The paper is 
organized as follows. 

In section 2, on the base of the Hopf-Cole solution, the asymptotic behaviour of the 
pulse at the nonlinear and the old-age stages is analyzed in general terms without assuming 
the detailed structure of the initial pulse. 

In section 3 the evolution of simple N-waves with regular and random amplitude is 
discussed. 

Section 4 deals with pulses with monochromatic carrier and section 5 with pulses with 
noise carrier. It is shown that the parameters of the asymptotic waveform depend weakly 
on the fine structure of the initial pulse, but that the old-age behaviour is very sensitive 
to the properties of the carrier. 

2 Solution of Burgers' equation and its large time 
asymptotics 

Our starting-point is the Burgers' equation 



which governs the propagation of plane acoustic waves in nonlinear dissipative media. 
Here v is the velocity of fluid and v is the fluid viscosity. Burgers' equation (|l|) has the 
Hopf-Cole solution [|16],[17j 



dv dv 



d 2 v 



(1) 



dt V dx 



dx 2 




(2) 



where U(x, t) is the solution for the linear diffusion equation 




(3) 



— oo 
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Here G(x,y,t) is the Green function of the linear diffusion equation 



G(x,y,t) 



\J Anut 



cxp 



and Uq(x) is the initial condition for this equation 



U(x,0) = U (x 
which corresponds to the initial condition 

v(x, 0) 



cxp 



(x - yf 

Aut 



2v 



(4) 



(5) 



v {x 



dipo (x) 
dx 



(6) 



The main goal of the present paper is the investigation of the old-age behaviour of 
the localized pulses. But first we discuss shortly the asymptotic evolution of nonlocalized 
waves. Consider a group of perturbation with the bounded initial potential |V>oOe)| < 00 
assuming that ipo(x) is a periodic signal with a period L Q or homogeneous noise with rather 
fast decreasing probability distribution of the potential ipQ. For such a perturbation in 
U(x,t) (H) we separate a constant component U: 



U(x,t) =U + U(x,t) . 



(7) 



Inserting (|7]) into (0) we see that U does not depend on time, since the y-integral of 
G(x,y,t) is unity. Here U(x,t) is a field with zero mean value (on the period or statisti- 
cally for noise). As times goes on, the viscous dissipation and oscillation (inhomogenity) 
smoothing causes the amplitude (variance) of the field U(x, t) to become less. At times 
when it amounts to \U\ <C U the solution (0) is equal to 

U x (x,t) 



v(x, t) = —2u- 



U 



(8) 



As U (x, t) satisfies the linear diffusion equation, then v(x, t) also at these times fulfils the 
linear equation. This testifies precisely to the fact that the evolution of the velocity field 
enters the linear stage. The accumulated nonlinear effects are described in this solution 
by the nonlinear integral relation between the initial velocity field Vq(x) and the fields 
U(x, 0), U (^]||), and are characterized by the value \Ax/jq\/u ~ Rcq. 
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Here Aipo is the characteristic change in amplitude of ipo, and R&o is the initial 
Reynolds number. 

From (H) it is easy to get the well known result, that for Re^ 3> 1 the initial harmonic 
waves asymptotically has also harmonics form, but with the amplitude not depending on 
the initial amplitude @.|3|]. 

At large initial Reynolds number the homogeneous Gaussian field vq(x) at the non- 
linear stage transforms into series of sawtooth waves with strong nongaussian statistical 
properties [2l,|25|. Nevertheless at very large time, when the relation (||) is valid, the dis- 



tribution of the random field v (x, t) with statistically homogeneous initial potential ipo{ x ) 
converges weakly to the distribution of the homogeneous Gaussian random field with zero 
mean value, and universal covariance function KHJ. But the amplitude of this function 



is nonlinearly related to the initial covariance function of the field vq(x) and increases 
proportionally exp(i?eg) with increasing of initial Reynolds number Rcq pT|,[23|. 



Let us consider the opposite situation when the initial perturbation is localized in 
some region. We assume also that initial potential ipo{x) vanishes very fast for x — > oo, 
so we have 

ip (x) = , \x\>L m , (9) 

where can be considered as the length of initial pulse. The condition (|9|) is identical 
to the assumption that the integral over the initial pulse is zero. 

It should be stressed that if the integral over the initial pulse is A?/> 7^ 0, then the 
initial perturbation transforms to the unipolar pulse with area A^o asymptotically. For 
Rcq = Aipo/v ^> 1 the form of this pulse is close to the triangular with the width of shock 



8(t) ~ ut/x s much smaller than the length of the pulse x s = y2\Aip \t. For such a pulse 
S(t)/x s ~ R&Q 1 = const and does not depend on time. Thus the nonlinear stage of the 
evolution prevails for all times. 

Let us now consider the time asymptotic behaviour of initially localized pulse which 
satisfies (0). If we rewrite @ as 

+00 

U(x, t) = l+ f [U (y) - l]G(x, y, t)dy , (10) 



then we can see that because of @ the integrand at right hand side of fllPf) is zero outside 
the region \y\ < L*. It is seen from (^) that the length scale ldif{t) of Green's function 
G(x,y,t) is 

W = • (11) 

If the length scale of Uo(y) is called Ljj, and the following condition 



l dif {t) > Lf/ 



(12) 



is valid, then Green's function can be considered as constant in the interval where the 
integrand at the right hand side of (|T0D is significantly different from zero. Thus we obtain 
from fllUp using (g) 

U(x,t) « 1 + G{x,y*{x,t),t)B . (13) 
Here y*(x,t), \y*\ < L#, is the value of y, in the neighborhood of which the integrand in 



' |10| ) gives the essential contribution in the integral. In (13) we introduce the notation 



B 



exp 



2u 



dy . 



(14) 



We assume B ^ 0. From (|^),(|l^),(|T4|) an approximate solution is obtained 

^^ Iff-t} . (15) 

t 1 + BG(x, y ,t) 

If the initial pulse is centered around y = 0, and the large scale of the field v(x, t) is 
much greater then L*, then we can put y* ~ and obtain 



. . x BG(x,0,t) 
v{x, t) = — - 



(16) 



t l + BG(x,0,t) 

instead of flTop. For large initial Reynolds number Rcq = A/2u, where A > is the 
maximum of ipo{y), the constant B flT4]) may be rewritten as the product of the maximum 
of the integrand and some length: L e ff 



B = L eff e A/2u , Re = A/2u 



(17) 
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Using the definition of Green's function (£|) and (|T7|) we can rewrite (|I~6| ) in the form 



exp 

X 


1 fx 2 A , 47TZ/t\" 

A + z/ln ^5— 




t 

1 + exp 




V 2 f L\„) 





In the limit of vanishing viscosity (z/ — > 0) we get from (|T8|) the well known solution for 
N-wave 



fx 



v(x,t) = < 
where x s is the position of the shock 



(19) 



j I oc | ^ j 



x s = V2M . 

At this stage the form and the energy of the pulse 

E(t) = [ v 2 (x,t)dx = 



(20) 



25/2^3/2 
3^/2 



; 2 i) 



are determined only by the value of maximum of initial potential A and does not depend 
on the fine structure of the pulse. This limiting case iy — > 0) for the pulses with complex 
inner structure was investigated in the paper |3T| . 

For large but finite Reynolds number the solution (|Tj|) , completed with a solution valid 
at x near x s and exhibiting the shock structure, is valid only at finite time. For finite t 
values and sufficiently small x we still have 



BG(x,0,t) > 1 . 



Then for 



\x\ < x q 



\ 



2tlA-u1n 



Anvt 



7-2 / ' 
L effJ 



(22) 

(23) 
(24) 



BG(x a ,0,t) = 1 

we obtain the increasing part (in x) of the N-wave solution (|19|). On the other hand, for 
x — > oo, t fixed, we find that v(x,t) fades away exponentionally. 
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For t — > oo, x arbitrary finite, we have 

BG(x,Q,t) < 1 



and the solution (|Lq ) can be approximated as 

x 



v(x,t) = B 



exp 



x 



Avt 



Using 



we can rewrite the condition of the validity of the solution fl26f) as 

A/2u 



VAthA > B 



\ 



At: ut 



T 2 



> e 1 



(25) 



(26) 



(27) 



Because ( pip gives a solution of the linear diffusion equation, the condition Q2"7| ) defines 
the linear stage of the pulse evolution. On this stage the pulse energy is 

j2 



£(t) = 5' 



7/ 



8vrf 3 



(2? 



Thus at the linear stage of evolution of the pulse has a universal form (|2q) and is 
determined only by the constant B, which is defined by the initial perturbation by the 



relation (14). The case B = can be understood by comparing the definition of B 
according to the equations (14) and (5) with the formulas (30) - (33) in ref. [19]. By 
this comparison it is clearthat the case B = corresponds to the absence of the Fourier 
component with n — 1 in the terminology of ref. [19]. The Fourier component with n = 
is absent here already for B ^ 0. The case B = excludes the N-wave solution (19). 

Below we will consider three cases of initial perturbation, assuming that the initial 
potential may be written in the form 



ip (x) = M(x)F(x) 



(29) 



where M(x) has the scale L — L*. First we consider the simplest case when F is 
constant, either deterministic or random value. Here we discuss the main properties of 
the wellknown solution (see ref. p0| ) for the N-wave. The cases of monochromatic and 
noise carrier will be considered in Sections 4 and 5. In these cases we assume that the 
scale Iq of the carrier F(x) satisfies the condition l Q L and then 



v(x,0) = v (x) ~ M(x)f(x) , f(x) 



-F'(x) 



(30) 
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3 The decay of a simple pulse with random initial 
amplitude 

3.1 The evolution of a regular N-wave 

We discuss the evolution of a one scale pulse, whose potential has the structure in (|29|): 

ipo(x) = M(x)A. 

Here A is a constant and M(x) is a dimensionless function with the scale Lq: 



(31) 



M(x) « 1 - 



x 

2Z§ 



+ ... 



(32) 



In particular we will consider the cases where the initial perturbation fl3"Tp, ( 52" ) is exact 
in the interval \x\ < V2L and ^(aO = outside this interval. We first assume A > 
and find from (Q) that the initial velocity pulse is an N-wave: 

, , xA 



V {X) 



L 



, Fl 



< V2L C 



v (x) = 0, \x\ > v2Lc 



(33) 



From (15D we then obtain 



fo(x) = exp 



Vo(^) 
2i/ 



exp 



- (i 



a; 



x < 



£/o(y) = 1, x > V2L C 



(34) 



From (|34[) we find that for this special case the length scale of Uo(y), first introduced 
in ([TT|), is (the square root of the inverted coefficient of x 2 in (|34])) 



,2i/ ,i L 
Lu ^ L (—)z = — r 

A Re| 



Re n = 



A 
2v' 



(35) 



The calculation of U (x, t), from which we obtain v(x, t) by use of (Q), can now be done 
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using (|10D with insertion of (|3~4 ) and 



[/(x,t) = 1 + 



V2L 



1 

1 + 2/ .,x 12 



y/Anvt J-V2L0 
1 



cxp 



1 + 51 





■/ 


< erf 






A 



T 2 I 



cxp 



A 1 



2z/ 

A A x 2 



2LI 



2v AvL\ M 
T 2 



1 exp 



{x - yf 

Avt 



dy 



V2L 



x 



1 + 



— < erf 
2 I 



V2L 



At 



x 



'Avt 



+ erf 



\Aut ( 1+ Ll 
V2L + x 
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+ erf 






A 



1 + 



At 
If/ 



J_f Af 



The definition of the error function is: 



(36) 



erf z 



7T JO 



exp(— t 2 )dt. 



The asymptotic behavior of the error function for large arguments is 



erf z = ± f 1 e 



r(|) r(§) 



+ ... 



z — > ±00. 



Using (H) we find from @§) for |x| < v^o, t finite, z//A < 1: 



U(x,t) 



'1 + 



At 



exp 



A 
27 



V 



x 



2^(1 + y^: 
^0 



(37) 



(3? 



(39) 



On the other hand, for \f2L$ <\x\ < \/2Lo 1 + -r^ ) we nave f° r * finite and v/A <C 1, 
vt/LiQ -C 1 and AtjL\ ^> z//A (the last inequality is necessary for the cancellation of the 
last two error functions in (36)): 



U(x,t) 



'1 + 



At 
T 2 



cxp 



A 
27 



V 



x 



2^(1 + ^: 



+ 1. 



(40) 
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/ L 2 \ 

For .'' • » 2A ( t + — ^ j , t finite and v j A < 1 we can neglect the second term at the 
righthand side of (|40D . Taking v — ► with t finite in (|39D or (|40| ) we obtain using (||): 

At 

3y < '\ «2/ e 



x L 2 
v(x, t) « — 



1 + 



At 



7~ 2 



f (x, £) = 0, |x| > x s 



(41) 



where x s is the shock coordinate 



x s 



\ 



2 -M t + i 



(42) 



However, for growing t-values the sharp discontinuity of the N-wave solution fliTD is 
broadened and we can calculate a shockwidth which depends on a small but finite value 
of v. Using ( |40D in (|2|) and assuming that 



(43) 



we obtain 



f (x, t) 



X 



exp 



A x 2 
2v ~ lui 



t L 



exp 



A x 2 
2u~ lui 



(44) 



+ 1 



Using ([14] ) and Q32|) we find for the present case in ([17]) L e jj = L J Attu / A and 



5 



/47rLgz/ 



A 



exp 



2z/ 



-^01 



' 2tt 
i?e 



3 _Re 



(45) 



and thus, using (f|), we see that the result (|44]) is a special case of (|16|). 

From the expression ( f4*4"|) we will now find the width of the shock. If we define the 
position x s of the shock as the coordinate where the wave amplitude has decreased to half 
its maximal value, we find from (|23| ) the following expression for x s : 



X s 



\ 



2t [A - v In 



AT 
LI 



(46) 
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Evaluating the expression ( fi4"D for x in the neighborhood of the shock, or more precisely, 
for 

\x - x s (t)\ 2 



vt 

we find 



< 1 (47) 



KM) = ^(l-tanh^-^V (48) 
where the shockwidth is given as 

In order to decide which of the two waveforms ([11]) and (f48|) is most appropriate we 

compare the expressions for x s fl4*2|) and (f46|). From these formulas we see directly that 

dx 

the shock velocity (v s = — -) decreases faster for growing t with nonzero viscosity The 

at 

zero viscosity expression (^) can be used as long as the correction to the N-wave shock 
position x = \JlAt is much greater in (ff^) than in fl4"6|), which means 

Ln v , ( At\ , t A 
_ ; 

^0 / f-M.ni 



77 >T In 72 , tln^ < -t M ,nZ, (50) 



where the nonlinear time tM,ni is defined as 

L 2 

tM,nl = -J- (51) 

In the notation tM,nh M stands for the kind of pulse and nl stands for "nonlinear". 

The energy of a pulse is defined as the integral over the pulse length of the square of 
the fluid velocity. We will calculate the energies of some of the pulses studied above. For 
the N-wave (^T[) we obtain 

„ 2 5 / 2 A 2 
E(t) = J v 2 ( Xj t)dx = — ^. (52) 

3 i _____ 

Thus for t ^> tM,ni and fl5"U[ ) still valid we see from (^) that E(t) behaves like A 2/^2 ( ]2T| ) 
and thus depends on the initial scale only in the next order of t, i.e. (^LqA 1 ^ 2 /t 3 ^ 2 ^j . 
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The energy of the wave under the condition t ^> tM,ni but (|50|) no longer valid is 
obtained from fl4~4]). We obtain 



+00 



E(t) 







1 



exp 



'A__x>_\ n2 
2v Aut J 



L 

1 H — p= exp 



At 



'A x 2 " 
2v Aut 



dx. 



After transformation of the integral in (^) we obtain 



E(t) 



2v l l 2 B 2 

Tit 3 / 2 



1/2 



2"S 5 



where .£> is given in ([45]). We need the formula ( [|3^], formula (2.3.13.27)) 



e x — z 



rrfx = r(a) [$(z,a,p + 2) - (p +!)$(*, a,p + 2)] 



where $ is given by the power series 



(53) 



(54) 



(55) 



$(z, = ^(p + n) a z n } \z\ < 1. 

n=0 



(56) 



Using (|55]), with a = 3/2 and p = 0, in ( |5l| ) an energy expression, valid for t ^> t M , n [ and 
(cf. (|50|) ) tin 3> —tM,nh is obtained: 



o 1/2 r>2 



7T 



^ 1 , x , B 3 , 
_ ,2 -$ ^,-,2 



(57) 



V 7 !^ 2' ' V y/^Ut T 

From (|57| ) the energy of the wave in the linear region is obtained by keeping only the 
power zero in the series for $: 



E(t) 



(5f 



A t^ 2 ' 

which is a special case of (2~8) with 5 given by ([4"5]). After introduction of the "linear" 



time tM,Un, the validity condition of ( 58f) is written: 

-^== < 1, t > t M im = tMn/exp(A/z/) = exp(A/z/)^ = t Mnl e 2Re ° . (59) 
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For the negative A the pulse decays much faster than for A > 0. For v — > at t > tuni 
the initial pulse transforms into so called S-wave 0,0, and the energy becomes 

o5/2 t3 

E(t) = '—f , (60) 

which is independent of the initial amplitude \A\ of the pulse. Of course the reason why 
the energy of the S-wave decreases faster with t than the energy of the N-wave is that the 
length of the N-wave increases with growing t in contrast to the S-wave. 

At large Reynolds number \A\/2u ^> 1 the constant B fll4|) which determines the 
linear stage of evolution is 

B = -2 3/2 L (61) 

and much smaller than for positive A with the same Reynolds number. The region of 
validity of linear regime is obtained from ( p5|) and (|6lD as 

t^L 2 /is . (62) 

Thus for the pulses with negative A the linear stage begins much earlier than for the 
pulses with positive A, for which the linear time tuun (|59D increases exponentionally with 
the Reynolds number Reo = Aj2v. 

Nevertheless, we can see on this simple example the main nontrivial properties of the 
old-age behavior of the pulses in the case of large Reynolds number. 

In the limit of vanishing viscosity {y — > 0) the energy of the pulse ( 2T|) ,( p^) and the 



shape of the pulse ([T9|),(|20|),(p]) do not depend on the length of L asymptotically for 



t — > oo. They are determined only by the maximum of initial potential ipo (in our case 
- of the area of triangular pulse of initial pulse A). Nevertheless, at the linear stage the 
energy depends not only of A, but is also proportional to L\ (|53D . We can remark that 
the energy of initial pulse is proportional to A 2 /Lq. 

The other property that we stress is that the amplitude B flPTD,(fi5D of the wave at 
the old-age stage depends exponentially on the amplitude A of the initial pulse. For 
small wave numbers k the Fourier component C(k,t) of the pulse depends linearly on 
k: C(k,t) ~ b(t)k. At the linear stage the slope b(t) of the Fourier component does not 
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depend on time: b(t) ~ B = const. At the nonlinear stage the growing of the slope until 
the very large linear time (|59|) leads to an exponentionally large value (in \jv) of b(t) at 
the linear stage (cf. (f45]) ). 

3.2 The decay of an N-wave with random amplitude 



Let us discuss the evolution of N-wave (|33f) with random initial amplitude. We assume 
that the amplitude A in ( [31]) has a gaussian probability distribution function with a zero 
mean value 

W(A) = -^= e - A2 ' 2tJ l . (63) 



From fl3"3p it is easy to see that mean initial field (v (x)) is zero, and relative fluctuation 



of the energy fl52|) at t = is 



SE(0) = < g2 <°» - <f °» 2 = 2 . (64) 

v ; (e(o)) 2 

In the previous section it was shown that there is strong difference between the decay 
of the pulse with positive and negative amplitude A. If we introduce E + (t) - the mean 
energy of pulse with positive amplitude A (N-wave), and E-(t) - the mean energy of 
pulse with negative amplitude A (S-wave), then at time 

t > t M ,ni ^ Ll/a^p (65) 

from (|52|) , ([BHD and ( |6^) we have 

E + (t)/E-(t) ^ {a^t/Llf/ 2 ~ (t/t MM f 2 » 1. (66) 

This fast decrease of pulses with negative A will lead to the generation of a field with 
positive mean value (v(x,t)) from an initial field Vq(x) with zero mean value (vq(x)) = 0. 

At time t ^> t M , n i ([55]) we can neglect the influence of the negative pulse on the mean 
field and use the expressions dl9[),(|20f) for the velocity. In different realizations we will 
have the N-wave with the same slope v' x = 1/t and random position of the shocks. 
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Let us introduce the cumulative probability of the random amplitude 

H 



Q H {H) = Prob{A <H)= J W{A)dA 



(67) 



Then from ©gOD 



v(x,t) = x/t with probability 1 — Qh(x 2 /2t) 
v(x,t) = with probability Qu(x 2 /2t) 



(68) 



From ( [68]) we obtain for the mean velocity (v(x,t)) and its variance a 2 (x,t) = ((v(x,t) — 
(v(x,t))) 2 ) the following expressions: 

„2\\ 



X 



X 



X 



a 2 (x,t) = -Q H - [1-Q H \- 



t 2 



2t 



2t 



(69) 
(70) 

(71) 
(72) 

where erf (z) is the error function fl37|) . 

At small x < (a^t) 1 / 2 the mean field (v(x,t)) is half the regular field with A ~ a^, 
due to the fact that only in that half of realizations, in which A > 0, we have a relatively 
slowly decaying N-wave. The mean field (v(x, t)) (|71~D has not a clear stressed shock front 
and has a very fast decaying tail at x ^> y/o^t: 

A 



For the Gaussian distribution of A fl63|) we have from (|67j), (|6 

,*,, t )=£(i-rfyy), 



(v{x,t)) ~ — ^= — exp 



7T x 



(73) 



It is easy to see that both mean field (v(x,t)) and variance a 2 (x,t) have self-similar 
behaviour 



(v(x,t)) 



**y 



X 



a 2 {x,t) 



(74) 
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where V, oy are given as: 



Kto)=yi-erf(^= 



4M = ^(l-er?(^ 



(75) 



Due to the self-similarity the relative integral fluctuation of the field 

+00 

I o%(x,t)dx 



5v v 



+00 

/ (v(x,t)) 2 dx 

—00 



3.356 



(76) 



does not depend on t and is of the order of unity. 

The situation radically changes on the linear stage when the field is described by the 
equation (|26|) . At this stage we have the self-similar evolution of the field 



where Vu n (y) is defined as 



v (x, t) 



Vun{y) 



^ ( -4= 



^7T 



cxp 



" 2 



(77) 



(78) 



and i? is a random amplitude fl45|) with nonzero mean value. For the gaussian distribution 
of F in (E9|) we have from (HT 



+00 



(B) = J (e M2 ^ Re °/ 2 - l) dx. 
—00 

Here we introduce the effective Reynolds number Re^ 



(79) 



Re = {F\x)) 1 l 2 /2u = a^/2 



v. 



(80) 



From (fr5|) we see that mean value (5) does not depend on the fine structure of the carrier 
f(x) and is always positive. The positive mean value (B) is a result of generation of mean 
field at the nonlinear stage. Here we consider the case of large Reynolds number, where 
we can use the asymptotic expression for B (pE5]). The n-th moment of B will be expressed 



through the probability distribution function (p.d.f.) of A (p3|) 

D nA/2u 



(81) 
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Using the steepest descent method for Reo > 1 we have from 



i) ^ e " as * 2 (82) 



Thus we have very fast growing of the momentum with increasing n. From (|77|) one can 
see that at the linear stage both mean field and variance are self-similar, and that the 
relative integral fluctuation 5v 2 does not depend on t. But in contrast to the nonlinear 
stage (|76|) , in the linear stage the relative integral fluctuation of the field 5v 2 is extremely 
large for large Reynolds number: 

5v 2 = {B2) - { 9 B)2 * V< (83) 
(B) 2 2 v ; 

We have similar behaviour also for the relative fluctuation of the energy 5E(t) 
At the nonlinear stage at a^t ^> L\ we have from ([H]), (|63|), fl64]) 



6E{t) = Y^k ~ 1 " 3 ' 376 • (84) 

Thus at the nonlinear stage 5E(t) does not depend on the initial scale L and the 
variance of the p.d.f. of amplitude A. 

At the linear stage the fluctuation of the energy is determined by the fourth and 
the second moment of B, and from ( |8"2"D we have also a strong enhancement of energy 
fluctuation ( j34[ ) 

8E(t) = ^e 4ite o. (85) 



4 The evolution of a pulse with monochromatic car- 
rier 

In this section we consider the evolution of a pulse with monochromatic carrier 

ipo(x) = M(x)A cos k x, v (x) ~ M(x)a sin k x. (86) 

Here A = ao/ko, and we assume that the length scale Lq of the modulation function M(x) 
is much greater then the period Iq of the carrier (Z = 2ir/k ). Below we consider two 
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large parameters Re® — A/2u and Lq/Iq. For different ratios of these parameters we have 
different regimes of wave evolution. 

The pure monochromatic signal (M(x) = 1) is characterized by the nonlinear time 
t s ,ni = stands for "signal") and the linear time t s ^ in = l/vk^ [^TJ. At t s>n i <C 

t <C t S) u n (Reo ^> 1) the monochromatic wave transforms into the sawtooth wave with 
the slope v' x = 1/t. The shock amplitude AV = lo/2t, as well as the energy of the wave 
E = /q/12^ 2 , does not depend on the initial amplitude. The shockwidth 

6 = 2u/AV = Avt/l Q (87) 

increases with time, and is, at t ~ t s ^ ini of the same order as a period. For t ^> t s t i in 
we have the linear regime of evolution, where the wave form is sinusoidal again with 
exponentionally decaying amplitude v(x,t) = 2vk§ exp[— ukot 2 ]sm k^x. 

For a pulse with monochromatic carrier a large-scale component vi(x,t) is generated. 
The interaction of the high frequency component with the large-scale wave vi(x,t) will 
change the evolution of the carrier. 



4.1 The nonlinear stage of evolution at large Reynolds number 



Below we give a short summary of pulse evolution at v — * based on the paper [[H]] , and 
discuss the influence of finite dissipation on the evolution of large scale and high frequency 
carrier. 

The nonlinear effect leads to the generation of the large-scale component vi(x,t), and 
at t 3> t Stn i this component obtains the stationary form 

vi(x,t) = -M'(x)A, (88) 

which is equal to the form of simple pulse vq{x) = M'(x)A with the same initial potential 
i/jq(x) = M(x) A. The evolution of the large-scale component is characterized by the 
nonlinear time tu,ni = Lq/A (|5l|) , and ( |38D is valid for t < tu.ni- At time 

t int = (89) 
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the energy of large-scale and small-scale components are of the same order. At t ^> t s>n i 
the evolution of the large-scale component is equal to the evolution of a simple pulse with 
ip (x) = M(x)A. 

For M(x) = (1 — x 2 /2Lq), \x\ < V2L the evolution of the large-scale component 



is described by the expressions fl42|), ([52]). At t ^> t S:n i, in the limit v — > 0, the 

evolution of the large-scale component is determined by only one parameter A = ao/fco of 
the initial perturbation. 

At t 3> t Stn i the amplitudes AV^ of the shocks of the small-scale components do not 
depend on the initial amplitude: AV = lo/2t, but due to the interaction with large-scale 
component they have nonzero velocity V r ~ —AM'(x r ) = vi(x r ) (|88D , where x r is the 
initial position of the shock. The distance between the shocks increases with time as 

l(t) = l (l + At/L 2 ). (90) 

The collision of the shocks of small-scale components moving with constant velocity of 

fine structure with shocks of the large-scale structure (021) leads to decrease of the number 

8L 2 



of shocks (see fig. 3 from |HJ]). The last collision occurs at time t = —r^-tM,nii and after 
this time a pure N-wave remains. 

At finite Reynolds number the width of the shocks of the large-scale component in- 
creases with time (|87D. The linear spreading of small structure leads to the increase of the 
distance between the shocks (|90]). Comparing ( p7|) with ([Xj) we find that, if the initial 
Reynolds number satisfies the condition 



4 Re » §, (91) 
o 



2u u ir 



then the nonlinear structure of shocks will be conserved. This is because the relative 
shock width is 



£(t) = v_ L 
l(t) A I 2 



4 T IT- ( 92 ) 
o 



In the opposite case, under the condition 

K Re < Ll/ll (93) 
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the nonlinear structure will dissipate before the nonlinear distortion of large-scale com- 
ponent begins. 

The evolution of the large-scale component at finite Reynolds number will be described 
by the same equations as the evolution of a simple wave fl4"7|) - (f49|). Only the position 
of the shocks x s (t) will be determined by some other equation (p3|), where L e ff depends 
on the fine structure of the initial pulse. This dependence leads to the sensitivity of the 
old-age behaviour on the fine structure of the initial pulse. 



4.2 Old-age linear stage evolution of pulse with monochromatic 
carrier 

The final linear stage of evolution is described by the equation fl26|) , where the constant 
B ([n]) is now given by 

+00 

B = J {exp [Re M (x)cos k x - 1] dx}. (94) 

—00 

At large Reynolds number Re^ = A/2v the constant B may be written in the form 
B = L e ff exp[i?eo] (0). In fact for the large Reynolds number the main contribution in 
the integral in ( |94"D comes from the neighborhood of points y r = l^r (Iq = 2ir/ko, r = 
0, ±1, ±2, ...). An evaluation of B by the steepest descent method then gives 

B = £ exp[Re M(y r )}. (95) 

r k jRe M(y r ) 



From (|95|) we see that the prefactor L e ff in (|17[) strongly depends on the ratio of two 
large numbers: the Reynolds number Re^ and the number (l = lnjk^). When the 

relation ([TI|) is valid, then only the first term r = in fl9"o] ) is significant, and from (|TTD, 
(|95|) we have 



Leff = -^d==- (96) 



k 

In the case of moderate Reynolds number, when the condition (93) is valid, we have 
from (0,|5D 

= (97) 
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The results above should be compared with expression ([45]) of B for the simple pulse 
for which we have L e ff = Lq^tc/ Reo. One can see that the modulation leads to faster 
transformation of the wave into the linear regime of evolution and consequently to decrease 
of the amplitude of the wave (B ~ L e ff exp[A/2u}) at the old-age stage. 

5 The evolution of a pulse with noise carrier 

In this section we will study the evolution of a pulse with noise carrier f(x) = —F'(x) 
([30]). We assume that the potential F(x) is homogeneous gaussian noise with rapidly 
decreasing covariance function 

B4p) = (F(x)F(x + p)) = a 2 R(p) = a 2 (l - ^ + ^ + . (98) 

In the limit of vanishing viscosity the continuous homogeneous field vo(x) = f(x) 
transforms into sequence of sawtooth pulses with equal slope v' x = 1/t and with random 
position of shocks. Due to the collision and merging of the shocks, their number decreases 
and the characteristic scale of random field increases. This effect makes all the statistical 
properties of the field self-similar and they are determined by only one scale l(t), which 
can be interpreted as characteristic distance between the zeroes of v(x,t) or between the 
shocks |lj,||,|2|. 

The evolution of the integral scale l(t) in time due to merging of the shocks ( see 
eq. 4.15, p. 170) 

l(t) = (M V2 ln- 1/4 (99) 
depends on only two integral characteristics of the initial homogeneous field 

°\ = (F 2 (x)), al = al/ll = (f 2 (x)). (100) 

Here Iq is the correlation length of the initial potential. Due to the merging of the shocks 
the energy density of the random wave (v 2 (x,t)) = l 2 (t)/t 2 ~ t^ 1 decreases slower then 
the energy of harmonic perturbation. 
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In case of the finite Reynolds number the thickness of the shock in the sawtooth wave 
originating from a monochromatic wave is given in (|8T|). We have the same expression 
for a random wave as well, where AV^ is the random amplitude of the shock. For the 
estimations we can assume AV ~ l(t)/t. 

The ratio between the integral scale l(t) and internal scale 5(t) is the effective Reynolds 
number Re(t): 

ft(1|a |,«M^ k -Wi), (101) 

0{t) V \tnnlJ 

where Reo and t n , n i ("n" stands for "noise") are defined as 

Re = a^/2u, t n>n i = Jq/o>. (102) 

Thus the effective Reynolds number Re(t) logarithmically slowly decreases with time, 
and the linear stage of evolution starts at very large time t n ^ in ~ t ntn i exp(Rel) ^> t ntTl i, 
(Re(t n ji n ) ~ 1). At this time the nonlinear effects become unimportant and the noise 
enters into the linear mode where its damping is mainly determined by linear dissipation. 
On the base of the solution (||), as it was shown in [^1|], the energy decays as 

al(t) = {vlljt^Re^e^l. (103) 

At the linear stage the distribution of the random field converges to the distribution 
of the homogeneous Gaussian field with zero mean velocity and variance according to 



01) 0. 

The evolution of the pulse with noise carrier in the limit of vanishing viscosity {y 



0) was considered in the paper | 31fl . It was shown that an initial perturbation v§{x) 
transforms to an N-wave. In the case, when the scale of the carrier Z is much smaller 
then scale of modulation function Lq, it was shown that the fluctuation of the shock 
positions is relatively small. 

Below we consider the properties of the pulse with noise carrier at the nonlinear stage 
for finite Reynolds number and the old-age behaviour of the pulse. 
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5.1 Nonlinear stage of evolution of pulse with noise carrier 

In the case of vanishing viscosity we can introduce two characteristic times: the nonlinear 
time for noise carrier evolution t njU i = Iq/ct^ and nonlinear time for modulation evolution 
tM,ni = L^/a^. At t ^> t n . n i the initial wave transforms into the sequence of sawtooth 
pulses, with the integral scale lM(x,t) ~ M 1 ^ 2 (x)lo(t/t n ,ni) 1 ^ 2 and the energy density 
E(x,t) = l 2 M (x,i)/t 2 ~ M(x)a 2 (t/t n ,ni)~ 1 depending slowly on the coordinate f31fl . At 
this stage the nonlinearity leads to partial depression of modulation and to generation 
of the mean field vi(x,t) = (v(x,t)) = —M'(x)a^(\nM(x)t/27rt nn i) 1 ' 2 pi] . Due to the 



merging of the shocks their number decrease, and at t 3> tu^ni the initial pulse with 
noise carrier transforms into an N-wave with random positions of zero and shocks. The 
position of the N-wave zero is localized in an narrow region A/ ~ Lo/(ln LI/2-kII) 1 ^ 2 near 



the center of the initial pulse |31| . The position of the shocks are determined by the 
equation (p0|), where A = ip m is the value of absolute maximum of the initial potential 
ipo(x) fl29l) . It was shown that the cumulated probability Qh(H) = Prob(-?/> m < H) has 



the form |31 



Q H (H)=exp[-N OQ (H)l 



H 



2 

(104) 



exp 



2*3 



Here N^ili) is the mean number of intersections of level H by the initial potential ipo{ x ) 
in the interval \x\ < (H), L is a characteristic scale of the modulation function M(x) 
([32]) (L ~ £*), = crv/ cr f i s an integral scale of the Gaussian homogeneous carrier 
potential F(x) and is its variance fllOOj ). 

The mean velocity (v(x,t)), with v(x,t) given by (|68|) and its variance a 2 (x,t) are 
described by equations (|69D , (|70|), where now the cumulative probability Qh(H) is deter- 
mined by equation ( |104| ), and not by the error function in (|7T|). It is easy to see from ( [74] ) 



that both mean field and variance have self-similar behaviour. From (|6^) , ((7(|), ( |104j) we 
have (cf. (74)) 

V(y, N) = y(l - q N (y)), a 2 v (y, N) = y 2 q N (y)(l - q N (y)), (105) 
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Qn(v) = Qh(v o>/ 2 ) = ex P 



V2N ( y 



r 



exp — — 



(106) 



V 2 V 4 

where y = xj y/cr^t, and iV is a large parameter, proportional to the number of correlation 
lengths on the whole extension of the pulse: 

N = -^U. (107) 



Contrary to the case of the simple initial pulse (|75|) , in the case of pulse with noise 
carrier the mean field and its variance have two different scales for N ^> 1 (see figs. 1 
and 2). The two scales in the "b" and "c" curves in the figures are the width of the shock 
and the length of the pulse. In the "a" curves only the length of the pulse remains. For 
N ^> 1 the function gjv(y) — at y < y* ~ (41nA^) 1//4 . The function qN(y) increases 
rather fast to 1 in the narrow region {y — y*)/y* — (41niV) . 

Thus for N ^> 1 the mean field has the N-wave similar structure with the dimensionless 
shock position 

/ T2 \i/4 

y. = y. = (4]nJ\0 1 / 4 =(2]n^|J . (108) 
The relative width of the shock of the mean field is 

^1 = *L ~ _J_ = I (109) 

x s y* ~ AlnN 21n{L 2 /27rl%)' 1 ; 

In the neighborhood of the shock position we can introduce a new variable z: 

y = y.(l + z/yt), (110) 

and from ( [106| ) we see that the shape of the front is described by double exponential 
distribution 

qn(v) = qd((y-y*)yi), q d (z) = e~ e \ (in) 

The variance a% is different from zero in a narrow region ( |109| ) near the shock position 
and for the relative integral fluctuation of the field ( [76] ) we have from ( |105| ), ( |106| ) 

s v i ~ ~ ^_ < i. (H2) 

y* 41niV v y 
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A finite Reynolds number leads to a finite width of the shock. The shock structure in each 
realization is described by the expressions fl4"8|) , fl4"5|) . The shifting of the shock position 
x s from the zero viscosity limit x s = ^2ip m t (A = ip m ) is given in (|23"D . At the nonlinear 
stage we can neglect this shifting. Then in each realization we have self-similar evolution 
of the pulse, and the relative width of the shock 5/x s = 2v/ip m does not depend on time. 
While for N ^> 1 (Lq ^> Zo) the maximum ip m of the initial potential is localized in the 
narrow region AH/H ~ QnN) near H ~ a,/, (2 In iV) 1/2 ( |T| formulas (107), (113)), we 
can estimate the relative width of shock front as 

x s - H ~ Re (lniV)V2' Ke ° ~ 2v i n6 ) 



Comparing ( |113| ) with (|109|) we see that the influence of finite viscosity on the mean field 



is unimportant if 



L 2 N 1/2 



ln 7r4 ^(lniV) 1 / 2 « J Re . (114) 



2tt/ 



In the opposite case, for extremely large ratio L /l , the width of the mean field will be 
determined by the viscosity. 

The displacements of shock position x s ( ^3|) from the vanishing viscosity position 



x s = y/2tp m t leads finally to the depressing of nonlinear effects. For t ^> tu n ~ L\^e^ m l v j v 
(cf.(P7|)) we have the linear stage of evolution. While for the random perturbation we have 
a large fluctuation of tu n , the cumulative action of nonlinearity, which is proportional to 
ti in , leads to strong fluctuation of the field at the linear stage. 

5.2 Old-age linear stage of evolution of a pulse with noise carrier 

At the old-age stage, when the evolution of the pulse is described by linear solution (|2~6D all 
the properties of the wave are determined by the constant B, given in flUp. The potential 
ijjo (x) (|29|) is a random Gaussian function. The mean value of B is given by the equation 
( [79|) , and does not depend on the fine structure of the carrier. For Re = a^/2u ^> 1 we 
have from ( |8"2"D 

(B)=L V2^-^e Re2 °/ 2 . (115) 
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Let us first consider the case of relatively large Reynolds number, when the integral 
(|I4]) may be calculated using the steepest descent method, and only the contribution of 
the absolute maximum is important: 



B= fm e ^ r (116) 

Here ip m = ipo( x m) is the value of absolute maximum of ip , an d = iPo( x m) is the 
second derivative of the potential at this point. To find the statistical properties of B we 
need to know the joint probability distribution of ip m and ip'^. From (|98|) it is easy to 
find the correlation coefficient for %p m and ip'^. r = If we consider the conditional 

probability distribution function Wfy'^/ipm), we then obtain for the conditional mean 
value (V'm)V'm an d variance J : 



While B increases exponentionally with ip m only the maximum ip m 3> cr^ gives a significant 
contribution to the average characteristics of B. Thus in (|116|) we can substitute (ip!^)^ m 
( |1171) instead of tp'^, and so we have 



B = l J-e^. (118) 

V i>rn 

From (|104j) we have for the probability distribution function of the absolute maximum 
of the initial potential ipo(x): 

W H {H) = 9 -§^ = -N'nWe-"-™ (119) 



This function is localized near 



and has a tail for H — > oo: 



H = a v ,(21nA^) 1/2 (120) 



Wh(H) = —e'^^l. (121) 
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For the calculation of the moments of B we can use the steepest descent method using 
the tail of the probability distribution function ( |121| ): 

(5 n ) = Lo/o- 1 f^) n/2 ^e^o/ 2 . (i 2 2) 



We have the main contribution for the n-th moment ( |122| ) at the point H* n ~ na^j2v = 
n Rcq. Comparing H* n with H given in ( |1 20|) we see that the inequality (|1 14|) is necessary 
for flTgp to be valid. 

Comparing ( |122| ) with the relation fl82| ) for the moment of B for a simple pulse we see 
that for the small-scale noise carrier the moments n > 2 depend on the integral scale Iq 
of the carrier. 

The relative integral fluctuation 5v 2 of the field (|76| ) on this stage is expressed for two 
first moments of B as 

h ^W-^rM 4 ' (123) 

and when the condition ( |114| ) is valid it is very large in contrast to the nonlinear stage 
( |112j ), where 5v 2 <C 1. For the noise carrier with the scale Iq <C Lq the relative integral 
fluctuation is the small factor lo/Lo times the fluctuation of the simple pulse with random 
Gaussian amplitude d83|). 

The calculation of (B), (AB 2 ) may be done directly on the base of the integral (|bj|), 
and for large Reynolds number we have the same equation (|123|) . For (1 /Lq) <C exp[— Re 2 ] 
we have relatively small fluctuation of 5v 2 , and with increasing of Iq/L the probability 
distribution of B approaches slowly the normal distribution (normalization). This normal- 



ization is similar to the normalization of the homogeneous field at the old-age stage f23 



But while the moment B n increases exponentionally with n, we have weak convergence 
to the normal distribution of B like in 1 12311 . 



6 Discussion and conclusion 



We have investigated the evolution of pulses with complex structure in nonlinear media, 
described by Burgers' equation. The investigation for vanishing viscosity was done in our 
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paper |3T| . There it was shown that the asymptotic form of an arbitrary initial pulse with 
zero area is an N-wave. It was also shown that the shock positions of the N-wave are 
determined by only one parameter of the initial perturbation - the value of the absolute 
maximum of the initial potential i/jq(x) (vq(x) = —dipo(x)/dx). It was also shown in the 



paper |JT| that for the noise carrier the fluctuations of the shock positions are relatively 
small if the carrier length scale Iq is much smaller than the modulation length scale Lo- 
in the present paper we have considered the evolution of pulses with complex structure 
at large but finite Reynolds number. On the base of the Hopf-Cole solution it is shown 
that the finite viscosity leads to a finite shock width 5, which does not depend on the 
fine structure of the initial pulse and is fully determined by the shock position in the zero 
viscosity limit. The other effect is the shift of the shock position from the zero viscosity 
limit position. This shift depends on the fine structure of the initial pulse, and as a 
consequence the time tu n , at which the nonlinear stage of evolution changes to the linear 
stage, is determined not only by the value of the maximum of the initial potential but 
also by the fine structure of the pulse. Because the amplitude of the pulse at the linear 
old-age stage is determined by the time tu n , the old-age amplitude is also sensitive to the 
inner structure of the pulse. 

In this paper the evolution of a simple N-pulse with regular and random initial ampli- 
tude and of pulses with monochromatic and noise carrier is considered. It is shown that 
the nonlinearity of the medium leads to the generation of a non-zero mean field from an 
initial random field with zero mean value. It is also shown that, at the nonlinear stage, 
the relative fluctuation of the field and its energy is of unit order for simple pulses and 
small for pulses with complex inner structure (Zo << Lo). However, at the old-age linear 
stage this fluctuation increases exponentially with increasing initial Reynolds number. 
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Figure 1: Self-similar mean velocity field form V(y) for the pulse with random initial 
amplitude ([751) (curve a) and for the pulse with noise carrier ( |105| ) for iV = 100 (curve b) 
and N = 1000 (curve c). 
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Figure 2: Self-similar fluctuation form <Jy(y) for the pulse with random initial amplitude 
( f75|) (curve a) and for the pulse with noise carrier ( |105| ) for iV = 100 (curve b) and 
N = 1000 (curve c). 
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